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Abstract 

Recently the CCFR Collaboration reported the measurement of the Gross-Llewellyn Smith 
sum rule Jq dFi^^'^'^^{x, = 3 GeV^) = 2.50±0.018(stat). Subsequently Kataev and Sidorov 
analyzed the Q^-dependence of this sum rule and pointed out a discrepancy between the 
results obtained via integration of the NLO fits to xF^i^XyQ"^) and the purely perturbative 
prediction. We suggest an explanation of this disagreement and show that the result of the 
CCFR measurement of the GLS sum rule integral can be determined from the knowledge of 
the value of A^ only. 

The Gross-Llewellyn Smith sum rule provides one of the important bridges between the QCD 
improved quark-parton model and the constituent quark model. Moreover, the quantity 

G{Q') ^ \ I' dx {f;;\x, q') + f!'{x, q')) (i) 

is one of the few physical quantities calculated so far up to the NNLO of perturbative QCD 
Q. In § the CCFR Collaboration reported the result 

G{Q^ = 3 GeV) = 2.50 ± 0.018(stat.) ± 0.078(syst.). (2) 

In p] Kataev and Sidorov analyzed this result and addressed an important question of 
the Q^-dependence of this quantity. They pointed out certain discrepancy between the 
predictions of perturbative QCD for this quantity and what they call experiment. In this 
note we explain the reason for their observation and argue that in a consistent treatment 
of the latter there is no place for such a discrepancy. We also comment on the procedure 
used in to evaluate (p. However, as the publication contains only a very sketchy and 
incomplete description of this procedure, we frequently consulted the unpublished thesis [Q 
for details. 

We start by recalling the familiar fact that in QCD the quantity G{Q^) can be evaluated 
in two different ways. First, it can be written as a purely perturbative expansion in some 

iln the following ^3(2;, Q^) stands for the average (f!^P{x, Q^) + F^P{x, Q^)) /2. 
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renormalization scheme RS (we use a = ag/TT as an expansion parameter instead of as) of 
the form 

G{Q^) = 3 (l - a(RS) (l + ri(RS)a(RS) + r2(RS)a2(RS) + ••■)) (3) 

While ri is known aheady for some time from |0, r2 has been calculated, for massless quarks 
and in MS RS, only very recently in The factor 3 on the r.h.s. of the above relation is 
of crucial importance. In the quark-parton model it corresponds to the fact that there are 
3 constituent (valence) quarks in the proton, but in QCD it has much deeper meaning and 
is based (see for a recent review of this subject): on the validity of 

• equal time commutation relations of currents made from the quark fields, 

• the asymptotic freedom. 

Both of these two properties are so fundamental to QCD that we cannot question the factor 3 
in (^) without abandoning the perturbative QCD itself. On the other hand were it possible to 
really measure (|I]), the result would be, as emphasized in the original papers |[, extremely 
important for the verification of the very foundations of QCD. 

The second theoretical route to G{Q^) exploits the nonsing let QCD fit to the depen- 
dence of xF^{x,Q'^), which, beside theoretical ingredients, requires also the specification of 
some initial condition. In this initial condition at some Qo "was taken in the form 

xF^{x, Ql) = Ax"(l - x)^ + Cx^ (4) 

while Kataev and Sidorov employed a simpler form, corresponding to C = in the above ex- 
pression. Fitting A, a, (3, C, 7, together with Aj^, to CCFR data one can compute F^i^x, Q^) 
at any x and and therefrom evaluate G{Q'^) at any The main observation in 0] con- 
cerns the disagreement between these two evaluations of G{Q'^). Similar observation has in 
fact been made already in 0] , but there the reason for the appearance of such a discrepancy 
turns out to be different than in 

In practice (|l]) is very difficult to really measure. Because of experimental limitations 
the accessible range of x varies with and prevents the measurement of G{Q'^) separately 
at each value of Q^. At low only very low x are accessible while for large only data 
at large x are available. Any experimental measurement of G{Q'^) therefore inevitably in- 
volves extrapolations to unmeasurable regions and must consequently be considered as a 
combination of data and some kind of theoretical prejudice. In p, ^ this extrapolation pro- 
cedure actually exploits the result of a nonsinglet QCD fit to F^{x, Q'^)\ So the experimental 
measurement of G{Q'^) actually assumes the validity of perturbative QCD. There is nothing 
illegal on this assumption, but, as emphasized above, we have to be consistent and include in 
the fitting procedure the constraint leading to the factor 3 in (^ as well. While the original 
determination p|, || does include this constraint on the QCD fit, the analysis of [§] it does 
not. 

Although the thesis Q contains many essential details on the determination of G{Q'^), it 
fails to provide convincing and full description of the NLO QCD fits used in the evaluation 
of G{Q^ = 3 GeV^) at so small Q^. First, there is the question of the number /i/ of massless 
quarks used in the evolution equations. As these are written for a fixed number of massless 
quarks, one must employ some kind of approximate procedure for crossing of quark mass 
thresholds. This is important in particular for the region of low Q^, below the charmed quark 
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threshold, which gives nonneghgible contribution to G{Q^) but which is certainly ineffective 
at = 3 GeV^. On the other hand some formulations in ||^ indicate that the fits used in 
determination of G{Q^ = 3 GeV^) used data with > 12 GeV^ only! In this case n/ = 4 
would be a good approximation, but then the fits do not actually fit the data at around 
3 GeV^! Although our argument is independent of the quality of these fits, this point should 
certainly be clarified. 

The problem of crossing the quark mass thresholds is only one of the subtle points in 
determining the GLS sum rule for so low Q^. There are others as well: 

• Target mass corrections. These can be easily included via the procedure of 0, but 
were ignored in extrapolation procedure of ^. 

• Higher twist terms. For = 3 GeV^ these turn out to be very important ||10|| ! 

• Cabbibo mixing. Below the charmed quark threshold the d — c is ineffective, thereby 
reducing the quantity G{Q^) by a factor cos^-i^c- 

In PI the reason for evaluating the quantity GIQ"^) at so low was motivated by the fact 
that low X region gives dominating contribution to it and 3 is just the evarage value of 
for the lowest x value measured in the CCFR experiment. However, looking at the data 
in this bin (see Fig. 4 below) it is clear that these data are obviously incompatible with 
the pure QCD evolution used in the extrapolation! Whether this is due to experimental 
problems, or physical effects not included in the pure perturbative QCD is another matter, 
but it is certainly much safer to follow the procedure employed in 0. There only > 12.6 
GeV^ was used in NLO QCD fits to CCFR data on ^3(2;, Q^), thereby avoiding most of the 
problems mentioned above, and G{Q'^) obtained by integration of these fits was compared 
to the purely perturbative expression (H) taken to the LO, with the couplant aiQ"^) defined 
to the NLO (for reasons discussed below) and with n/ = 4. The main observation of 
concerns the discrepancy, displayed in Fig.l, between these two results and in particular its 
Q^-dependence. In the rest of this note we identify the source of this discrepancy and show 
how to avoid it. 

Recall that as a consequence of the factorization theorem in QCD the non-singlet (NS) 
structure function F^i^x^Q"^) can be written as a convolution 

F,{x, Q') = C %^s{y. M)G^s (jT, -, a{f,)] (5) 
Jo y \M y J 

between the NS distribution function gNs(x,M), satisfying the evolution equation 

dgNs(a;,M) dy 



dlnM /„ >ste)P.s ^.a(M)) (6) 



and the "hard scattering cross-section" G{Q/M, z,a{fj.j). The scales M, /i appearing in (^) 
are the so called factorization and renormalization scales and the couplant a(M) satisfies 
the familiar equation 

-ba^ (1 + ca(M) + ■•■), (7) 



dlnM 
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where the first two coefficients b,c are unique functions of nf. The branching function Pns 
as well as the hard scattering cross-section Cns admit the following perturbative expansions 



P^s{z,a{M)) = a{M)P{,'i{z)+a\M)P^^^[z) + 
C^s{Q/M,z,a{fi)) = 6{z)+a{fi)CP{Q/M,z) + --- 



>{i)/ 



(8) 
(9) 



In the NLO approximation only the ffist two terms in each of the expansions (0) ^ (H) ^ (H) 
are taken into account. Turning the above expressions into the corresponding ones for the 
moments defined as 

F,{N,Q')= /'x^F3(x,Q2)dx (10) 
Jo 



and analogously for other quantities, we get an explicit expression for F3{N,Q'^) ITT 

FsiN, Q2) = gNs(iV, M)Cns(Q/M, N, a{fi)) 



where 



q^s{N,M) = AN 



a(M) 



1 + ca{M) 



[1 + ca{M))-''^^' 



l + a(/.)|4^1n^ + ^ 



:iii 



:i2) 



(13) 



and d^^ , i = 0,1 are the moments of the LO and NLO branching functions P^^^' 



P 



NS ' 



z) and 

respectively, and the hard scattering scale fi is in general different from the factor- 



NS 



ization scale M. 



Recall also that while d^^^ is unique, d^P 



''N 



is essentially arbitrary, defining 
the so called factorization convention (FC) |Tl|, 0. Two invariants appear in (pAl): the 



quantity kn, which is independent of both M and d^j^\ but depends for general on the 



choice of the renormalization scheme of the couplant |T^; and the constants A^, which are 
independent of anything. As argued by Politzer |Tl|], these constants must be kept fixed 



when the factorization scale M or factorization convention, defined at the NLO by d^^\ are 
varied! Moreover, these constant provide alternative way of specifying the necessary bound- 
ary condition on the solution of the evolution equation. Note, however, that fixing An does 
not correspond to specifying the initial condition on F3(x, Q^) at any finite Qq. Rather the 
constants A^ determine the asymptotic behaviour of F^lN, M^) as M — oo. 

For = 0, corresponding to the sum rule (|T|), the situation is particularly simple. 



as do'^ = and in the FC employed in both p and (called universal) also c/q^'' = 



(1) 



Moreover the value of Kq is unique and equal to —1. Consequently, q^si^ 
Q^-independent and (|TT|) reduces to 

Fs{0,Q') = GiQ') = Aoil - a{fi)) 



A 



IS 



(14) 



which coincides with the purely perturbative expression (^ in the LO provided Ao = 3 and 
the same couplant a(/i) is used in both formulae. In the general case, when /i ^ M, the 
couplant a(/i) used in (|]) can be taken only to the LO approximation, but if we identify, as 
is the usual practice, jj, = M and use the same couplant in both (||), (|^), a(M) in (|T4D is in 
the NLO approximation. 
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The crucial point of this note is the observation that for the quantity G(Q^) the basic 
principles of QCD imply Aq = 3! Taking this into account in the NLO fits to the data on 
Fs{x, Q^), the integrals over these fits must coincide with the results of purely perturbative 
predictions (^! There is simply no place for the kind of discrepancy observed in 0. This 
conclusion holds independently of the quality of the NLO QCD fit to the data. Note also 
that while in the universal FC the constraint = 3 is equivalent to the condition 



dxgNs(a;,<5^ 



(15) 



for all Q^, in general FC this is not the case and (|T^) may depend on Q^. 

The discrepancy observed in is thus a direct consequence of the fact that in their NLO 
fits Kataev and Sidorov did not impose the constraint Aq = 3. The technique employed in 
to write down the QCD predictions for the evolution of F3{x,Q^) is based on the 
idea 



T3| , |Tj, |l5l of expanding the convolution (g) in the set of Jacobi polynomials 6 



a/3. 



X 



orthogonal on the interval (0,1) with the weight x"(l — x) 



F3{x,Q'' 



X 



XI 



fc=0 



a/3. 



xia 



al3 



(16) 



where the "Jacobi moments'' 
mentsF3(j,Q2) 



af3 



(Q^) are given as linear combinations of conventional mo- 



a/3 



(17) 



j=0 



In the boundary condition, necessary to specify such a convolution, is taken in the form of 
the initial condition (|) with C = 0. From the fit of ^3(0;, Q^) to the CCFR data at all > 
12 GeV^ its authors then determine the values of the parameters y4,a,/?, Aj^, describing 
F3{x,Ql), and evaluate G{Q^) at Q^. Repeating this procedure with Ql corresponding to 
each of the 14 values of at which the data of exist, yields the points with error bars in 
Fig. 1. Their comparison with the purely perturbative LO results (continuous curves) reveals 
the discrepancy mentioned above. In principle even a single fit is sufficient to evaluate G{Q^) 
at all Q^, but to integrate the expansion (0) multiplied by 1/x down to x = is potentially 
dangerous because of the oscillatory character of the Jacobi polynomials. To integrate the 
initial condition (^ is, on the other hand, a straightforward and safe matter. It must, 
however, be kept in mind that athough all these 14 fits to F^^XjQ'^) are equally legal NLO 
approximations, they lead in general to numerically different results at all Q^. 

To prove our claim we have repeated the analysis of 0, evaluating for each initial Qq the 
value of the constant Aq. We have used the same method of Jacobi polynomials, described 
above, as in . Technical aspects of our calculations (number ot terms taken into account in 
the expansion (|r^), accuracy achieved etc.) can be found in |T^. For the class of conditions 
(I) Aq can easily be expressed as a function of A, C, a, /5, 7 



An 



A 



r(a)r(/3 + i) , c 

T(a + p + l) 



+ 



1 



1 



l-a(Qo) ' ' 

Qo in (HID- Turned around the condition = 3 implies a constraint between 



settmg /i 

the mentioned parameters, for instance 



A 



And 



am) - - 
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r(Q + /3 + i) 

T{a)T{l3 + l) 



(19) 
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In Fig. 2 we plot r{Ql) = Ao{Qo)/3, as given in (113), with A,a,P determined from fits of 



(p!6|) to CCFR data, using the parametrization (H) with C = 0, (i.e. the one used in 0), as 
a function of Qq. We have repeated this exercise also for three indicated cuts on used 
in the fits. For all of them the same pattern of the dependence Ao{Qq) emerges: for low Qq 
Aq is below 1, while for large Qq it rises above 1. Note that ^(Qo) is equal to the ratio of 
the two results in Fig. 1. For all chosen samples there is thus a value of Ql for which the 
unconstrained fit is consistent with the condition Aq = 3. This point depends somewhat on 
the cut, but lies roughly around 100 GeV^. The corresponding values of and x^/d-f. 
are plotted as functions of Ql in Fig.3. 

In the next step we imposed (for the simplest case C = 0) on our fits the constraint 
(110), guaranteeing thus Aq = 3. The dashed curves in Fig.3 describe the results of these 



constrained fits. From Figs. 2,3 we conclude that 

1. the value of extracted is only weakly dependent on the initial Qq, but 



2. depends sensitively on the Q^ 



cut 

3. for unconstrained fits x^/d.f. is an mcreasmg' function of Qq, which 

4. decreases with increasing Q^ut- 

5. As expected, the constrained fits yield essentially the same values of Arg^ and /d.f. 
for those Qg which r{Q'^) = 1. For lower, as well as higher, values of Qq the fits 

MS- 



have substantially higher x^/d-f- and yield dramatically different A|^. 



The last point demonstrates the importance of a proper choice of the initial Qq if an uncon- 
strained fit is performed and its correlation with the assumed form of the initial condition (^) 
and the value of Q^nt parameter. To further illustrate the difference between the constrained 
and unconstrained fits and indicate the importance of the choice of Qq, we compare in Fig. 
4 their results to CCFR data for the sample with Ql^^ = 12 GeV^. The values of extracted 
the values of the extracted A^ (see also Fig. 3) between 260 — 270 MeV are in reasonable 
agreement with the result Aj^ = 237 ± 36 MeV, obtained in 0]. 

In principle we could use more general form of the initial condition than in (^), adding 
further parameters and looking for such forms which would exhibit flatter dependence r(Qg) 
than in Fig. 2. In our own investigations we have taken slightly different form of the 
additional term (to the conventional Ax°'{l — x)^) in the initial condition (^) and namely 

Fs{x, Ql) = Ax'^il - xf{l + 7x) (20) 

The corresponding results are not significantly different and we therefore show only some 
typical examples in Fig. 3a. The basic feature of these results is, as expected, lower x^/d-f- 
of the fits. More detailed investigation of this question certainly makes sense but would take 
us beyond the immediate subject of this note. 

Summarizing the preceding considerations we claim that provided 



the NS fits are performed with the constraint Aq = 3, i.e. for instance ([19D 



• and the LO purely perturbative approximation to (|^) uses the same NLO couplant 
a(M) as the NLO fits to ^3(3;, Q^) 
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these two ways of evaluating the GLS sum rule G{Q'^) must give exactly the same result. 
The discrepancy noted in is entirely due to the fact that in their fits Kataev and Sidorov 
do not impose the constraint Aq = 3. 

Finally, let us briefly return to the original determination of the GLS sum rule in 
which contains basically the same observation as in 0. There is, however, a subtle difference 
between the procedures employed in these two papers. In |^ the value G{Q^ = 3 GeV^) = 
2.50 ± 0.018 is obtained by extrapolating the NLO fits to = 3 GeV^, fitting this ex- 
trapolation to the form (^) with C = and integrating the resulting analytical formula. 
Using the LO term of (j^), but with LO couplant a(M) only, they get, on the other hand, 
2.65 ±0.03. They speculate that higher order terms and/or higher twist effects are, perhaps, 
responsible for this marked difference. In fact the reason is much simpler and again related 
to the constraint Aq = 3. As already emphasized, in the FC used in this constraint is 
equivalent to the condition 

f\Ns{x,Q^)dx = 3 (21) 
Jo 

which is explicitly imposed on their fits. Had they integrated the extrapolated form of their 
fits at = 3 GeV^, they would get exactly the result (0) with a^^^{Q), the NLO couplant. 
Had they furthermore used the purely perturbative expression (j^) to LO, but with a^^'~'{Q), 
their two results would coincide. However, by fitting the extrapolated form of -^3(0;, Q^) to (^) 
(with C = 0) there is no more any guarantee that the resulting fit will satisfy the constraint 
Aq = 3. And indeed, it does not\ Substituting the values of A, a, /3, published in Q to 
(p!8|) we find Aq = 2.73! This, together with the fact that the purely perturbative result 
with the NLO couplant equals (for the fitted value of Ay^ = 237 ± 36 MeV) 2.71 ± 0.03, 
explains quantitatively most of the effect ((2.73/3) x 2.71 = 2.46). Note that the value 
^Ms ~ 210 ± 28 MeV mentioned in [Q, comes from fit to a combined F2 (at large x) and 
data. 
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Figure captions 



Fig. 1: The quantity G{Q'^) as obtained by purely perturbative expansion (^) (solid 
and dashed lines) compared with its evaluation via integration of non-singlet fits to 
F^[x, Q^) (points with error bars). Dashed lines correspond to la error of Aj^. Taken 
from the ref. 0. 

Fig. 2: Qq dependence of the quantity r{Ql) = Aq/S extracted from nonsinglet fits to 
FsiXjQ"^) for three different cuts on and two options on 7 in (^). Typical errors 
(not shown) of r{Ql) are 0.03. 

Fig. 3: and x^/d-f- corresponding to nonsinglet fits with variable initial Ql for the 
same three classes of fits as in Fig. 2 and 7 = 0. In a) also the results corresponding 
to variable 7 are displayed. Errors (not displayed) of are typically 40 MeV. 

Fig. 4: The results of our fits to CCFR data on ^3(0;, Q^) for two values of initial Qq = 
10, 100 GeV^, with or without the constraint Aq = 3. Only data at > 12 GeV^ are 
included in the fits. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9410292v2 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9410292v2 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9410292v2 
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